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Dynamics  of  Homogeneous- Heterogeneous  Reactors 


R.  J.  Olsen,  W.  R.  Williams,  X.  Song,  L.  D.  Schmidt  and  R.  Aris 

Department  of  Chemical  Engineering  and  Materials  Science 
University  of  Minnesota 
Minneapolis,  MN  55455 

Abstract — In  many  industrial  processes  both  heterogeneous  and  homogeneous  reactions  can 
occur,  and  these  reactions  are  strongly  coupled  through  heat  and  mass  transfer  and  reactive 
intermediates.  Surface  reactions  can  heat  the  boundary  layer  sufficiently  to  ignite  homogeneous 
reaction,  can  deplete  the  boundary  layer  of  a  limiting  reactant,  and  can  act  as  a  source  or  sink 
of  free  radicals. 

We  have  studied  combustion  reactions  over  catalytic  foils  and  in  catalytic  monoliths  for 
simple  oxidations  such  as  CH^  +  O2,  NH3  +  O2,  CaHg  +  Oj,  CH<  +  NO  and  NH3  +  NO. 
For  most  reactions  and  geometries,  we  observe  heterogeneous  and  homogeneous  ignitions  and 
extinctions  and  self-sustaining  autothermal  steady  states.  The  range  of  compositions  for  which 
an  autothermal  state  exists  depends  sensitively  on  the  inlet  gas  velocity. 

For  stagnation  point  flow,  reaction  over  a  catalytic  foil  can  be  modeled  by  a  one-dimensional 
boundary  value  problem  in  which  saddle- node  bifurcations  correspond  to  ignitions  and  extinc¬ 
tions.  Computations  for  a  model  with  simplifled  kinetics  reproduce  the  experimental  dependence 
on  input  power  to  the  foil  and  inlet  gas  composition.  Calculations  also  show  that  the  aynamics 
are  affected  strongly  by  the  inlet  gas  velocity. 

1.  Introduction 

Many  processes  of  commercial  importance,  such  as  catalytic  oxidation,  involve  simultaneous 
homogeneous  and  heterogeneous  reaction  pathways  (Satterfield,  1980;  Pfeffeile  and  Pfefferle, 
1987).  Examples  of  catalytic  oxidations  include  the  production  of  syngas  from  methane,  the  <ur 
oxidation  of  NH3  over  platinum  to  produce  NO,  and  the  Andrussov.  process  for  the  synthesis  of 
HCN. 

The  nature  of  the  coupling  between  the  homogeneous  and  h'*terogeneous  pathways  must 
be  understood  if  homogeneous- heterogeneous  reactors  are  to  be  used  effectively.  The  coupling 
occurs  through  heat  and  mass  transfer  between  the  catalyst  and  its  surroundings  and  through 
reactive  intermediates  that  are  common  to  both  pathways.  The  surface  reaction  can  heat  the 
boundary  layer  sufficiently  to  ignite  homogeneous  reaction,  can  deplete  the  boundary  layer  of  a 
limiting  reactant,  and  can  act  as  a  source  or  sink  of  free  radicals.  While  surface  reactions  are 
frequently  observed  to  exhibit  ignition  and  extinction,  such  phenomena  have  received  less  at¬ 
tention  in  coupled  systems.  The  determination  of  ignition  and  extinction  points  is  an  important 
step  in  characterizing  the  dynamical  behavior  of  a  reactor  since  these  points  provide  information 
about  operating  limits  of  the  process. 

The  catalytic  combustion  of  premixed  CsHs/air  mixtures  in  stagnation  point  flow  over  plat¬ 
inum,  which  will  be  considered  here,  is  a  typical  homogeneous-heterogeneous  reaction  system. 
The  kinetic  parameters  for  global,  one-step  rate  expressions  of  Arrhenius  form  have  been  deter¬ 
mined  for  both  the  homogeneous  and  the  heterogeneous  reactions  (Westbrook  and  Dryer,  1981; 


2505 


2506 


R.  J  Olsen  et  al. 


Dll 


Hiam  et  ai,  1968).  Using  these  rate  expressions  along  with  known  values  for  the  physical  prop¬ 
erties,  we  are  able  to  compare  experimental  and  calculated  bifurcation  behavior  as  a  function  of 
fuel/air  ratio  and  power  input  to  the  catalyst.  Based  on  the  agreement  between  our  calculated 
results  and  experimental  data,  we  have  extended  our  previous  analysis  (Song  et  ai,  1991)  to 
study  the  effect  of  additional  parameters  on  the  behavior. 

The  momentum,  energy  and  mass  balance  equations  for  plane  stagnation  point  flow  admit  a 
change  of  variables  that  aUows  for  the  model  equations  to  be  written  in  terms  of  a  single  spatial 
coordinate.  For  steady  state  conditions,  the  model  equation  is  a  one-dimensional  boundary 
value  problem  with  2n  -f-  5  variables  where  n  is  the  number  of  chemical  species.  The  use  of 
global,  one-step  rate  expressions  retains  aU  aspects  of  the  coupling  except  for  the  role  of  radical 
intermediates. 

2.  Dynamical  Behavior  in  Propane  Oxidation  over  Pt 

Schematic  one-parameter  bifurcation  diagrams  at  four  different  feedstream  compositions  are 
shown  in  Fig.  la.  Each  curve  represents  a  typical  experiment:  the  inlet  gas  temperature,  mass 
flow  rate,  and  inlet  gas  composition  are  fixed  and  the  input  power  to  the  catalytic  foil  is  varied, 
using  an  experimental  apparatus  shown  in  the  inset  of  Fig.  2a.  The  measured  response  is  the 
temperature  of  the  surface  at  the  stagnation  point  (Ts).  At  low  and  high  power,  curves  are 
essentially  linear  due  to  simple  resistive  heating  of  the  foil.  For  intermediate  values  of  power, 
the  response  is  nonlinear  due  to  the  exothermicity  of  the  reaction,  which  has  a  self-accelerating 
effect  on  the  rate.  In  this  intermediate  region  a  transition  from  a  kinetic-limited  reaction  to  a 
mass  transfer-limitM  reaction  occurs. 


Power  ^  %  Fuel  ^ 

Figure  1:  Schematic  representations  of  one-parameter  (power)  and  two-parameter  (power  and 
%fuel)  bifurcation  diagrams.  Each  curve  of  (a)  is  obtained  by  varying  power  with  all  remaining 
parameters  fixed.  A  family  of  curves,  of  which  four  are  shown,  is  generated  by  varying  %fuel. 
1  he  surface  temperature  at  the  bifurcation  and  autothermal  points  for  each  one-parameter  curve 
is  plotted  as  a  function  of  %fuel  to  produce  the  two-parameter  bifurcation  diagram  of  (b). 

At  the  lowest  fuel/air  ratio,  Ts  is  a  single-valued  function  of  power  (as  illustrated  by  curve 
A)  because  the  heat  liberated  by  the  reaction  does  not  influence  the  system  significantly.  At  a 
critical  fuel/air  ratio,  the  bifurcation  diagram  has  a  vertical  tangent  as  shown  in  curve  B.  Further 
increases  in  the  fuel/air  ratio  lead  to  multiple  steady  state  behavior.  An  ignited  state  that 
sustains  itself  in  the  absence  of  input  power  is  called  an  autotherm.  The  onset  of  autothermal 
behavior  is  shown  in  curve  C,  where  the  zero  power  line  is  tangent  to  the  extinction  point. 
Another  feature  shown  by  curve  C  is  homogeneous  ignition,  indicated  by  the  vertical  arrow. 
Experimentally,  homogeneous  ignition  typically  results  in  a  flame  which  flashes  back  to  the 
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reactor  inlet,  although  stable  flames  in  the  reacting  boundary  layer  have  been  observed  (Williams 
et  al,  1991).  Bifurcation  diagram  D  exhibits  a  stable  autothermal  state  and  terminates  in  a 
homogeneous  ignition  point. 

Ignition,  extinction  and  autothermal  behaviors  are  sketched  as  a  function  of  inlet  gas  com¬ 
position  in  Fig.  lb.  This  type  of  diagram  is  the  result  of  projecting  the  ignition,  extinction 
and  autothermal  points  from  each  of  the  individual  curves  in  Fig.  la  on  the  (%fuel,Ts)- plane. 
We  have  chosen  this  representation  of  the  two-parameter  behavior  instead  of  the  more  familiar 
projection  onto  the  (power, %fuel)-plane  in  order  to  emphasize  the  dependence  of  the  surface 
temperature  on  the  fuel/air  ratio.  The  boundary  of  the  oval  region  in  Fig.  lb  represents  the 
ignition  and  extinction  points  for  the  heterogeneous  reaction.  The  interior  of  the  oval  region  cor¬ 
responds  to  the  unstable  steady  states.  The  piecewise  linear  representation  of  the  autothermal 
surface  temperature  as  a  function  of  composition  follows  from  the  assumption  that  the  kinet¬ 
ics  in  this  state  are  limited  by  mass  transfer.  If  a  single  stoichiometry  describes  the  reaction, 
the  autothermal  temperature  will  attain  iis  maximum  value  at  the  stoichiometric  fuel/ air  ratio. 
The  final  feature  of  Fig.  lb  is  the  curve  of  homogeneous  ignition  points,  found  at  relatively  high 
surface  temperatures  and  over  a  fairly  narrow  range  of  composition. 

3.  Experiments  and  Calculations  for  Propane  Oxidation 

The  experimental  apparatus  consists  of  an  electrically  heated  platinum  foil  suspended  in  a 
Pyrex  tube,  as  shown  in  the  inset  of  Fig.  2a,  through  which  the  reactant  gases  flow  upwards.  The 
lower  section  of  the  tube  is  packed  with  glass  beads  to  insure  that  the  gas  stream  is  well-mixed 
and  is  wrapped  with  heating  tape  so  that  the  inlet  gas  temperature  may  be  varied.  Surface 
temperature  is  measured  either  by  an  optical  pyrometer  or  a  Chromel-Alumel  thermocouple 
spot-welded  to  the  surface.  The  ignition  and  extinction  points  for  the  surface  reaction  are 
determined  by  locating  the  limits  of  steady  state  multiplicity  on  the  surface  temperature  vs. 
power  curves.  Homogeneous  ignition  is  defined  by  the  appearance  of  a  flame.  Ignition  and 
extinction  points  are  generally  reproducible  to  within  ±  20 °C  and  are  not  very  sensitive  to  the 
orientation  of  the  foil  with  respect  to  the  direction  of  gas  flow. 

The  experimentally  determined  surface  temperatures  for  heterogeneous  ignition  and  extinc¬ 
tion,  autothermal  behavior  and  homogeneous  ignition  are  displayed  in  Fig.  2a.  The  abcissa  R 
is  defined  as  R  =  $/($+  1),  where  $  is  the  fuel/air  ratio  divided  by  the  fuel/air  ratio  for 
the  stoichiometry  of  complete  combustion  to  CO2  and  HjO,  which  is  4%  for  propane  in  air. 
Except  for  fuel-lean  compositions  near  the  hysteresis  point,  the  surface  temperature  at  which 
heterogeneous  ignition  occurs  is  «200°C,  essentially  independent  of  composition.  Heteroge¬ 
neous  extinctions  take  place  at  400-500  °C.  These  occur  over  a  restricted  range  of  compositions 
due  to  the  predominance  of  autothermal  behavior  in  this  system.  To  a  good  approximation,  the 
autotherms  are  linear,  indicating  that  the  kinetics  are  mass-transfer  limited.  The  discontinuity 
in  the  slope  of  the  autotherm  in  the  fuel-rich  regime  coincides  with  the  visual  observation  of 
coke  formation  on  the  foil.  Homogeneous  ignition  is  observed  for  fuel/air  mixtures  between  2% 
and  8%,  with  a  minimum  surface  temperature  for  ignition  of  about  1000 °C  at  w3%. 

The  equations  for  planar  stagnation  point  flow  (see  the  inset  of  Fig.  2b),  written  in  dimen¬ 
sionless  form,  are 


(la) 

(lb) 

^+ySc^-Scr=0, 

(Ic) 
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Figure  2:  Homogeneous- heterogeneous  ignition  and  extinction  curves  for  CsHg  oxidation  over 
Pt.  The  experimental  data  are  shown  in  (a)  and  the  calculated  curves  in  (b).  The  experimental 
apparatus  is  sketched  in  the  inset  of  (a),  and  the  geometry  used  to  derive  the  model  equations 
is  shown  in  the  inset  of  (b).  The  model  accurately  reproduces  the  compositions  corresponding 
to  the  onset  of  multiple  steady  state  behavior,  autothermal  behavior  and  homogeneous  ignition. 


The  results  obtained  for  the  model  are  given  in  Fig.  2b.  The  qualitative  agreement  with  the 
experimental  data  is  remarkably  good  in  excess  Oj  where  the  rate  parameters  for  the  heteroge¬ 
neous  reaction  are  reliable.  The  surface  temperature  for  heterogeneous  ignition  is  reproduced 
accurately,  as  are  the  inlet  gas  compositions  for  the  onset  of  steady  state  multiplicity,  autother¬ 
mal  behavior  and  homogeneous  ignition.  The  maximum  autothermal  temperature  and  homo¬ 
geneous  ignition  temperature  occur  at  the  stoichiometric  composition  for  complete  combustion, 
reflecting  the  underlying  global  rate  expressions  incorporated  in  the  model. 

4.  The  Effect  of  Inlet  Gas  Velocity 

Given  the  success  of  the  model  in  reproducing  the  system  behavior  as  a  function  of  input 
power  and  composition,  we  continued  our  study  of  the  model  behavior  by  examining  the  effect 
of  varying  other  parameters.  We  describe  here  the  bifurcation  behavior  with  respect  to  inlet 
gas  velocity,  retaining  input  power  as  the  primary  bifurcation  parameter.  Kr,  hi,  Pw  and  the 
dimensionless  reaction  rates  r  and  rg  all  depend  on  the  inlet  gas  velocity. 
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Figure  3:  The  effect  of  inlet  velocity  on  the  homogeneous-heterogeneous  ignition  and  extinction 
of  CaHg  oxidation  over  Pt  for  a  composition  of  3%  propane  in  air  and  a  reference  velocity  vq  of 
5  cm/s.  Curves  in  (a)  and  (b)  show  Ts  as  a  function  of  normalized  power  (Pw  =  (v/vo)Pw)  at 
the  indicated  velocities.  The  curves  in  (c)  show  the  dependence  of  Ts  on  velocity.  Homogeneous 
ignition  occurs  in  the  absence  of  heterogeneous  ignition  at  sufficiently  low  velocities.  As  the 
velocity  is  increased,  the  homogeneous  ignition  is  suppressed,  and  at  very  high  velocities,  the 
autothermal  state  also  disappears. 

Fixing  the  composition  at  3%  propane  in  air  and  the  inlet  gas  temperature  at  300  K,  we 
obtain  the  bifurcation  diagrams  of  Fig.  3.  The  quantity  v/vq  is  the  ratio  of  the  inlet  velocity  to  a 
reference  velocity  of  5  cm/s,  which  was  the  velocity  used  in  Fig.  2b.  In  Fig.  3c  the  bifurcation  sets 
and  autothermal  curves  have  been  projected  on  the  (Ts,  v/vo)-plane.  Six  different  bifurcation 
diagrams  are  observed  as  the  inlet  velocity  is  varied.  Of  these,  the  diagram  at  very  low  velocity 
exhibits  no  bifurcations  and  the  diagram  at  very  high  velocity  exhibits  a  simple  hysteresis  loop 
resulting  from  heterogeneous  ignition  and  extinction  only.  Bifurcation  diagrams  at  intermediate 
velocities  are  displayed  in  Figs.  3a  and  3b.  For  the  lowest  velocity  shown  in  Fig.  3a,  there 
is  homogeneous  ignition  in  the  absence  of  heterogeneous  ignition.  For  a  velocity  one  order  of 
magnitude  larger,  a  narrow  range  of  hysteresis  arising  from  the  heterogeneous  reaction  is  visible 
along  with  the  homogeneous  ignition.  Homogeneous  ignition  without  heterogeneous  ignition 
and  autothermal  behavior  was  not  found  in  the  variation  of  inlet  composition  discussed  above. 
The  bifurcation  diagram  of  Fig.  3b  for  v  =  Vo  is  common  to  Fig.  2b  (R  =  3/7).  Increasing  the 
velocity  one  order  of  magnitude  further,  we  find  that  homogeneous  ignition  is  suppressed  while 
heterogeneous  ignition  and  an  autotherm  persist.  The  range  of  velocities  shown  in  Fig.  3c  cannot 
be  realized  in  any  single  experimental  apparatus.  AdditionaUy,  effects  due  to  free  convection 
and  flow  in  the  sonic  regime  are  not  included  in  Eqs.  1-3.  However,  it  is  the  trend  in  qualitative 
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behavior  shown  in  Figs.  3a  and  3b  which  is  of  primary  interest,  and  the  model  is  reasonably 
accurate  for  this  range  of  velocities. 

5.  Conclusions  and  Summary 

We  have  shown  that  a  model  for  homogeneous-heterogeneous  combustion  of  CaHs  based  on 
simplified  rate  expressions  reproduces  many  of  the  interesting  features  of  the  observed  exper¬ 
imental  dynamics  as  inlet  composition  is  varied.  Having  verified  the  qualitative  accuracy  of 
the  model,  we  are  able  to  simulate  the  dynamics  in  regions  of  parameter  space  that  are  diffi¬ 
cult  to  probe  experimentally.  For  example,  by  varying  inlet  gas  velocity  over  several  orders  of 
magnitude,  the  sensitivity  of  the  model  to  this  psirameter  is  obtained.  Both  the  experimental 
and  computational  results  illustrate  that  the  presence  of  the  catalyst  has  a  marked  effect  on 
homogeneous  ignition  behavior. 

This  work  was  supported  in  part  by  the  Minnesota  Supercomputer  Institute  and  the  Army 
Research  Office  (contract  number  DAAL03-89-C-00S8  with  the  University  of  Minnesota  Army 
High  Performance  Computing  Research  Center). 

Notation 

All  quantities  listed  are  dimensionless 

/  =  stream  function 
6  -  temperature 
(3  =  heat  of  reaction 
mi  =  fuel  mass  fraction 
rs  =  heterogeneous  reaction  rate 
Kr  =  radiative  heat  transfer  coefficient 
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Tj  —  distance 

Pr  =  Prandtl  number 

r  =  homogeneous  reaction  rate 

Sc  =  Schmidt  number 

Pw  =  input  power 

hi  =  conductive  heat  transfer  coefficient 
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